Noncollinear magnetic ordering in small Chromium Clusters 
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We investigate noncollinear effects in antiferromagnetically coupled clusters using the general, 
rotationally invariant form of local spin-density theory. The coupling to the electronic degrees of 
freedom is treated with relativistic non-local pseudopotentials and the ionic structure is optimized 
by Monte-Carlo techniques. We find that small chromium clusters (N < 13) strongly favor non- 
collinear configurations of their local magnetic moments due to frustration. This effect is associated 
with a significantly lower total magnetization of the noncollinear ground states, ameliorating the 
disagreement between Stern-Gerlach measurements and previous collinear calculations for Cri2 and 
Cri3. Our results further suggest that the trend to noncollinear configurations might be a feature 
common to most antiferromagnetic clusters. 
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E£ ■ I. INTRODUCTION 

P 

I ' Magnetic,-properties of transition metal clusters have become the subject of intensive research, both from the exper- 
imentaLrlj and theoretical point of viewETH. One of the most interesting and challenging aspects of that field is the 
.£3 ' subtle interplay between geometric structure and magnetic ordering which has mostly been investigated for ferromag- 
netic 3d-clusters and 4d-clusters. Finite size effects and a reduced dimensionality rftep lead to a significantly different 

c/2 ■ magnetic behavior from the bulk with clusters having enhanced atomic momentaaO, larger anisotropy energieald or 
an altered temperature dependence of the magnetizationlij. 

Almost all theory to date employed the local spin-density approximation (LSDA) with the assumption that the spin 
density matrix is diagonal in some particular frame. In that special case the spins are automatically collinear along 
a fixed quantization axis. The only generalized spin-density calculation for clusters that treats the electron spin as a 
vector observable and a function of position has recently been performed by Car and coworkers 113. They have shown 
that noncollinear configurations exist in Fe3 and Fes although the effect on structure and energetics of these ferromag- 
netic (FE) clusters is not very pronounced. On the other hand an unconstrained orientation of the quantization axis 
is knowit-to play a key role in describing various non-ferromagnetic systems like the 7-phase of bulk ironEJ, disordered 
systemsE-3 or ultra thin Fe films_with a partial antiferromagnetic (AF) couplingcj. Furthermore, the work on Fe/Cr 
(refill) and Ag/Cr multilayersE£l demonstrated how the competition between AF ordering and frustration of the Cr 
moments leads to noncollinear arrangements in form of a spin-density wave. Although the importance ojLfrustration 
in AF systems seems to be evident and-was discussed in detail for embedded Cr clusters by Pastor et al.tj and more 
OS \ generally by Manninen and coworkers!^, the possibility of noncollinear effects has not yet been considered for AF 
clusters. 

c/3 \ In this paper, we present the first general local spin-density calculation for clusters of AF materials. Besides the 

» 1 motivation given above, we have decided to explore noncollinear effects in chromium clusters for two reasons: First, 

chromium is particularly challenging amongst the 3d-elements due to its maximal spin alignment in the atom (3d 5 4s 1 ): 

£^» all valence electrons have parallel spins, i.e. the total magnetization is S=3. This leads to a variety aLppique effects 
including an unusually short dimer bonding length or a repulsion for FE coupling at small distancesEJEj. The other 
reason why we chose chromium is the possibility to compare with the thorough ab-initio study of Cr clusters by 
Cheng and Wangpj. These authors recently employed the conventional LSDA for N at < 15, finding all clusters to be 
antiferromagnetically coupled. Their structures can serve as a benchmark to check our computations concerning the 

3 j ■ restriction to collinear configurations - an aspect that is crucial to extract the importance of noncollinear effects. 

In Section O we outline the fundamentals of our theoretical approach and briefly motivate the structure optimization. 
This involves the rotationally invariant LSDA to minimize the electrons and a relativistic, non-local pseudopotential 



for their interaction with the ions. In Section III we discuss some numerical tests and present our results. We 
find noncollinear spin-configurations for all investigated clusters. We demonstrate the influence of noncollinearity 
on the geometry and on the total magnetic moment and discuss how this reduces [the discrepancy concerning the 
magnetisation of some chromium clusters between the experiment by Bloomfield et alfi-3 and some previous theoretical 
resultaH3'E3. 



II. THEORETICAL AND NUMERICAL BACKGROUND 

The density functional theory in its most general form, as developed by Barth and HedinEj, allows the orientation of 
each spin to vary with position. The wavehmctions are described by complex two-component spinors \& = (\&o-, ^a- 1 ), 
where a and a' denote the spin indices, and the degrees of freedom are the elements of the single-particle spin-density 
matrix 

Paa'(r1=£^ )CT (f)** CT ,(f) . (1) 

i 

Assuming this matrix to be diagonal, the usual local spin-density functionals are parameterized in terms of p\i(r) =: 
p-\(r) and Pn(r) =: p\,(r) only In that special case, the spins are necessarily collinear along the chosen quantization 
axis and the exchange-correlation potential is obtained via V£ c = dE xc [pj, pj jdp a . However, rotational invariance 
requires that the true variables are the eigenvalues 7if(f) and n\(f) of the spin density-matrix p aa >{f). We can thus 
apply standard local spin-density functionals (we chose the formula of Perdew and Wangcll) by evaluating the potential 
in a locally diagonal frame. The transformation is carried out following the work of Kiibler et alB3 who used the 
spin-1/2 rotation matrix 

/ cos ^fl e i4>(f) s i n 2ifi e ^m \ 

U V) — \ . 9(r) Ism ofr) =±<h(r) ( 2 ) 

\ -sin -j 1 e 2 9{ ' cos -^ e 2 i"s> I 
to locally diagonalize the spin-density matrix: 

^ U aa (r) p aa , (r) U*, p (r) = S a/3 n a (r) (3) 

By working in this representation we express dE xc jdp aa i by E xc / dn a plus the introduction of local spin rotation 
angles <fr (f) and 9 (fV which are the local azimuthal and polar angles of the magnetization density vector. They are 
computed from eqn.(|3j) through the requirement of vanishing off-diagonal elements as 



(Re m (rO) 2 + (Imp u (r)) 2 



1/2 



^( r -)=-tan- 1 ^4g ) e(r)=t^^- 1— ^ (4) 

These new degrees of freedom complicate the mean-field equations and lead to an exchange-correlation potential V xc 
in form of a complex matrix in spin-space 

v xc = i (f t xc + Vf c ) i + i (y T xc - v^ c ) ff ■ d , (5) 

where d is a position dependent unit vector along the direction of the vector (Re/3||(r), Im. p^(r), Pii(f) — Pair))- 
The presence of the second term in the exchange-correlation potential allows a general coupling of the up- and down- 
compouejxts of the spinor wavefunctions. The wavefunctions belong to the class TSDW of the Fukutome classification 
schemecl To interpret the magnetic properties, we compute the vector magnetization density rh{r) by expressing the 
spin-density matrix in the form 

/3(f) = 0.5 n(r) I + m(r) ■ a . (6) 

We associate magnetic moments with individual atoms by integrating each component of m(r) within a sphere cen- 
tered on the ions, giving us the local magnetic moment vectors /2 a t. The integration radius is chosen to be one half of 
the smallest interatomic distance in each cluster to avoid overlap and the resulting spheres contain about 80 — 90 % 
of the magnetization density. Taking a larger radius for more distant atoms has a negligible effect on the orientation 
of their magnetic moment, although its magnitude becomes somewhat bigger. 

As usual, we only treat the valence electrons explicitly, taking care of th£L.ionic core with a pseudopotential approx- 
imation. We use the relativistic pseudopotential from Goedecker et al.EI, which contains a local part plus a sum 
of separable gaussians, optimized to represent a transferable non-local pseudopotential on a coordinate mesh. The 
multiplication of the wavefunctions with the non-local part can be limited to a small region around the ions as the 
radial projectors fall off rather quickly. However, 20 integrations within the covalent radius of each atom need to be 
performed to correctly account for non-local effects in chromium. The energetics at small ionic separations inside the 



clusters further requires us to include the 3s and 3p semi-core electrons into the variational space. Our pseudopotential 
also includes spin-orbit terms which fix the orientation of the total magnetization M to the ionic structure thus giving 
rise to magnetic anisotropy. The implementation of the L ■ S operator is not too costly because we have to deal with a 
complex spinor-structure anyway. Spin-orbit effects enable us to test the validity of the usually applied atomic sphere 
approximation by studying the intra-atomic dispersiontJ. 

We have carried out an unconstrained structural search by fully minimizing electronic and ionic degrees of freedom. 
To find the ground state and stable isomers, the ionic positions were optimized by Monte-Carlo sampling applying 
the technique of simulated annealing. After some Metropolis steps, the electronic wavefunctions are updated with 
Kohn-Sham iterations. The optimization of the ionic geometry involves a minimization of the one-ion energies and is 
explained in detail in refEx The static Kohn-Sham equations are-solved in a combined coordinate and momentum 
space representation by using an efficient damped gradient iteratiorH. Local operators are applied on coordinate space 
wavefunctions while the kinetic energy and the action of the spin-orbit operator are copputed in momentum space 
applying fast Fourier techniques. The Poisson equation is solved via the FALR-methocO. As it is more convenient 
for most physical observables, electronic wavefunctions and densities are stored on a three-dimensional coordinate 
space mesh. We perform our calculations in a cubic box with a mesh spacing of 0.32 a.u. and up to 64 grid points 
in each direction. We checked that the mesh siza-was big enough to avoid artifacts from the boundaries. A detailed 
description of our numerics can be found in refsEj. 

III. RESULTS AND DISCUSSION 

Before discussing our results we mention some of the various tests we performed in order to increase our confidence in 
the Hamiltonian and its numerical implementation. The dimcr plays a key role in the description of small chromium 
clusters. It is known that ite-ipwbtle electronic properties demand a high accuracy of the Cr-Cr interactions and 
the numerical representationc3c3. By applying the pseudopotential in the semi-core version, our binding energy 
(1.98 eV) and bonding length (do = 3.25 a.u.) for the antiferromagnetic ground state were in good agreement with 
experimental results^ (1.56 eV ± 0.3 eV, d OX p = 3.19 a.u.) and previous all-electron or pseudopotential calculations^! 
The correct energetic order of the single particle levels as a function of the intra-molecular distance, the symmetry 
of the wavefunctions and the properties of the ferromagnetic coupling (d = 5.2 a.u.) could be reproduced as well. 
The same holds for the bonding length of the CrO molecule which deviated from the experimental result by 1.3 %. 
Additionally, we achieved degenerate d-states up to a level of 1 % and the correct energetic order of 3d and 4s-levels 
in the chromium atom. The rotationally invariant spin-density theory was checked by letting the FE configuration 
of the dimer relax to the AF ground state. Like in the collincar theory, all spins were initially restricted to point in 
the z-direction. The wavefunctions and energies of the final result turned out to be identical with the ground state as 
computed in a separate collinear approach although the quantization axis of both atoms had rotated by ± 90 degrees 
during the iteration. This confirms the degeneracy of the electronic properties with, respect to the orientation of their 
spin. Furthermore we have been able to verify the result of Car and coworkerstJ concerning the noncollinear spin 
arrangement of Fe3 . 

Our results for the energetically lowest noncollinear configurations of small chromium clusters Ci'n (3 < N < 
6) are presented in Fig.fjl. All structures except the one of Cr4 represent ground states. The corresponding total 
magnetization is shown in Table |. Our geometric and magnetic structures are obtained by performing up to 50 full 
Monte-Carlo runs per cluster starting from arbitrary ionic coordinates. A couple of thousand Kohn-Sham iterations 
are usually necessary to completely relax the electronic degrees of freedom. This is because the numerical convergence 
with respect to the direction of the local moments (governed by competing interatomic exchange interactions) is much 
slower than with respect to their magnitude which is determined by stronger intra-atomic interactions. 
The principle effect which leads to noncollinear arrangements in chromium clusters can best be demonstrated in 
&3. Our calculation restricted to collinear spins gives a trimcr that basically consists of a dimer plus a loosely 
attached third atom, very similar to the result of Cheng and WangEa. It is obviously impossible for the atoms to 
couple antiferromagnetically with all their neighbors at a time so that a FE interaction between two atoms becomes 
inevitable - a situation that we refer to in the following as frustration. The repulsive effect of the FE coupling at 
smaller distances pushes one atom away and results in a large interatomic distance of 4.97 a.u.. If the noncollinear 
channel is accessible, however, the trimer is free to achieve the closest to AF-like coupling it can by rotating two local 
magnetic moments in the x-z plane (see Figiy). As the frustration is now reduced, the third atom comes much closer (d 
= 3.67 a.u.) and enables a buildup of molecular orbits which enhances the total binding inside the trimer significantly. 
The associated reduction of |/2 a t | for the formerly isolated atom results in a decrease of M at from 2 /zb to 0.69 /ib- The 
energy difference with respect to the collinear ground state AE nc is 0.083 eV per atom (see Table |) which amounts 
to 7.8 % of the binding energy of &3. A look a the interatomic distances and the tilted local moments suggests that 



the dimer lost its dominant role. However, the fact that the ground state is not an equilateral triangle with angles of 
120° between the local magnetic moments (this would equal the best possible AF-like coupling) indicates that some 
trace of the strong dimer binding from the collinear calculation still persists. Unlike the situation in Fe3, the linear 
isomer of Cr3 did not favor noncollinear spins. 

The noncollinear spin structure of Cr3 can be understood as a compromise between the energetically very favorable 
AF coupling and frustration. This can be seen easily with the very simple Hamiltonian H — P'Yl li< j <J% ■ (fj for three 
spins on an equilateral triangle. Here (3 is negative for AF coupling. In a collinear restriction, the lowest energy state 
formed by a product wavefunction is (TTI) with an energy expectation (H) = j3. Taking instead the state with 120° 
angles between the spin directions gives a lower energy of (H) = 3/2 j3. For higher atomic spins (as they occur in our 
numerical computations) the preference of the noncollinear configuration would of course become more pronounced 
due to the larger number of exchange interactions. 

In Fig.0 we show the rotation angle 9(r) in the x-z plane of the trimer including its contour lines. The azimuthal 
angle <p(r) is zero at any gridpoint due to the magnetic anisotropy energy although this quantity is rather small. But 
it seems that even a difference of a few meV per atom is enough to keep the magnetization inside the plane of the 
trimer, an observation which agrees with the findings of rcfO. The orientation of the magnetization density vector 
is remarkably uniform in the regions surrounding the atoms (9 ~ ±74°, 0° with a small oscillation of A9 — ±10°) 
and varies significantly only right between them where the charge density is essentially zero. This is why a smaller 
integration radius only influences the magnitude of /2 a t and not its orientation. The rapid change from 9 = —90° to 
6 = 90° is related to a spin flip in these interatomic regions. We find an intra-atomic spin dispersion of around 6° 
which partly comes from spin-orbit coupling. But the dispersion is also induced by the trend to an AF coupling with 
the neighbors. The same especially holds for the variation of the spin direction close to the domains of other atoms. 
The change of 9 at the ionic positions is an artifact of the pseudopotential approximation. A detailed analysis further 
shows that these features are common to all investigated chromium clusters. 

The collinear ground state of Cr4 was found to have a rectangular structure (not shown) with bonding lengths of 3.35 
a.u. and 4.62 a.u.. The double dimer structure of this geometry appears to be resistant to non-collinear effects as the 
distance of two atoms with parallel moments is very large (5.7 a.u.). But in the case of the lowest isomer, a rhombus, 
frustration becomes important again resulting in noncollinear spins (see Fig.pl) . The argumentation follows the case of 
Cr3, the triangles in Cr 4 are very similar and have only slightly larger bonding lengths. A higher total magnetization 
of M at = 1-33 fiB leads to an energy gain of 0.12 eV per atom with respect to the corresponding collinear state of the 
rhombus and halves the energy difference to the rectangular ground state. The isomer of Cr 4 is unique in that the 
energetically favored noncollinear configuration has a bigger total magnetization than its collinear counterpart. The 
magnitudes of /I a t, however, were almost identical in both cases. It is also noteworthy that the rhombus structure 
provides angles between the local moments of almost exactly 120° although the bondings are not equivalent. This 
seems to be related with the larger moments of the furthermost atoms (A|/2 at | = 0.6 /ib)- 

The properties of the noncollinear ground state of Crs can be understood if one considers the pentamer as consisting 
of three triangles. Although the bonding lengths are bigger, each of these triangles exhibits the same basic features as 
Cr3 and the local moments are again arranged in such a way that the best possible AF-like coupling is achieved. We 
also see that the magnitude of /2 at decreases with increasing coordination number. Note that the reflection symmetry 
of the geometry is the same as the symmetry of the magnetization. The gain in binding energy with respect to the 
best collinear state is 0.054 eV per atom (— 3.7% of the binding energy) which is less than for the trimer. But the 
total magnetization is again clearly reduced from 1.03 to 0.53 /^b per atom. The lowest isomer is a bipyramid with a 
noncollinear spin structure as well and M at = 0.79 /Ub- Its energy difference of 0.6 eV to the ground state is 0.25 eV 
higher than in the collinear case. CY5 is the only cluster we found where the noncollinear and collinear geometries 
differ by more than just a variation of the bonding lengths. Our lowest collinear state looks similar to the geometry 
of Cheng and Wang and has a 2v-symmetry. However, its total magnetization of 1.05 /ib is somewhat higher than 
their result of 0.93 ^b- 

The shape of Crg, on the other hand, resembles the collinear one very closely. Ctq is the smallest cluster that displays 
a fully three dimensional geometric and magnetic structure in the ground state. It consists of three dimers distributed 
over two triangles in which frustration sets in. Each of the /I a t is exactly antiparallel to the moment of its partner 
atom in the dimer. The azimuthal angles are cf> = ±19° for the atoms in the foreground and 4> = 0° for the ones 
in the background which are slightly closer. The bonding lengths in the triangles are about 6% shorter than in the 
collinear case and the dimer distances are somewhat bigger. It is important to note that in spite of a vanishing 
total magnetization and although the dimer seems to recover a certain influence, tilted spins are still energetically 
favored. However, the gain of AE nc = 0.022 eV per atom only accounts for 1.4% of the binding energy of Cr6- 
This indicates that the impact of noncollinear effects on the energetics becomes gradually smaller with rising cluster 
size. Further calculations for Cr7 and Crg confirm this trend and predict that the noncollinear geometries essentially 
equal the collinear ones for N at > 6. On the other hand our results so far show that noncollinear spin configurations 
considerably reduce the total magnetization of the ground states although it might not be reflected too much in the 



binding energy. This can be understood in terms of the competing interatomic exchange interactions which involve 
only small energy changes during the rotation of local moments (see above). 

We shall finally discuss how noncollinear effects might provide a way .-to ameliorate the disagreement between the 
measured total magnetic moments of Ctn (N > 9) by Bloomfield et al.Ea and previous theoretical resultsEIEl The 
Stern-Gerlach experiment extracted an upper bound for M at , M at < 0.77 hb, assuming a superparamagnetic behavior 
of the chromium clusters. Earlier calculations, however, report values for the magnetization of some clusters that are 
much higheita. The values from Cheng and Wang for Cri2 (M at = 1.67 (ib) and &13 (M at = 1.06 (ib) also exceed 
the experimental limit. The trend from smaller clusters gives one the hope that noncollinear effects might reduce the 
difference. However, an unconstrained simultaneous calculation of electronic and ionic degrees of freedom surpasses 
our computational resources for clusters as large as Cri2. Instead, we start the optimization procedure from the 
geometries of Cheng and Wang. This is a reasonable approach in the light of the very good agreement with our 
collinear structures and the small changes in geometry that are induced by noncollinear spins. Our final magnetic 
and geometric configuration of Cri2 is shown in Fig.pl The free relaxation leads to a shortening of the bonding 
lengths between the corner atoms of about 5 % and slightly bigger distances of the capping atoms but the bulk-like 
bec-structure of the collinear geometry clearly persists. The z-components of the local magnetic moments vary on 
alternating x-y planes, but only the moments of the corner atoms have significant x- and y-components. All the spins 
of the lightest-colored corner atoms point inside the cube towards the central atom, whereas the spins of the slightly 
darker corner atoms point outside and away from the next atom. This dispersion is related to some frustration of the 
corner atoms in connection with the preferred AF coupling to the moments of their nearest neighbors (darker atoms). 
A magnetic arrangement like that can be seen as a precursor to the bulk behavior in form of an antiparallel order 
between neighboring atomic layers. The tilted spins of the corner atoms (4> = ±43°, 9 = ±34° and (j> = ±45°, — ±30°) 
have pairwise opposite x- and y-components so that the sum of their net-moments in positive z-direction is reduced. 
The moments of the darker atoms, however, show almost no dispersion in negative z-direction. All this results in a 
much smaller total magnetization of M at = 0.81 /ib which is now very close to the experimental limit. The remarkable 
reduction of AM at = — 0.86/^b is associated with AE nc = 0.011 eV per atom which equals only 0.5 % of the binding- 
energy of Cri2. A very similar situation leads to M a t = 0.59 ^b in the case of Cri3, a value that is even below the 
experimental limit. The noncollinear gain amounts to AE nc = 0.008 eV per atom. The Ci'13 result is preliminary and 
might slightly change (AM at < ± 0.05 /ib) because the electronic degrees of freedom were not completely relaxed. But 
we can still conclude that the total magnetization of Cri2 and &13 as obtained with the general LSDA represents a 
considerable improvement with respect to the experiment. 

IV. CONCLUSIONS 

We present the first study of noncollinear effects in antiferromagnetically coupled clusters by applying the general, 
rotationally invariant LSDA for the electronic degrees of freedom. Their interaction with the ions is described in terms 
of a relativistic, non-local pseudopotential which has been thoroughly tested. The magnetic and geometric structures 
are obtained by employing a simulated annealing technique for the ionic optimization together with interlaced Kohn- 
Sham iterations which update the wavefunctions simultaneously. 

Our collinear configurations, which we use in order to extract the effect of noncollinear spins, agree very well with 
the results of Cheng and Wanga. We find that all investigated chromium clusters show a pronounced trend to 
noncollinear spin configurations. This is caused by a subtle interplay between the preferred magnetic order and 
frustration, a situation which can in principle occur in all clusters of elements that favor antifcrromagnetic spins. 
Therefore our considerations appear to be of a more general nature although the special properties of chromium 
indicate that noncollinear effects could be less dramatic in other transition metals. The results for N at < 13 show that 
the influence of noncollinearity on various observables becomes gradually smaller with rising cluster size. Induced 
changes of the cluster geometry are generally restricted to an alteration of the bonding lengths, with the exception of 
Crs. However, we observe a universal reduction of the total magnetization which is significant even in those clusters 
for which a variation of the magnetic distribution is not clearly reflected in the binding energy any more (Cri2 and 
Cri3). This effect is related to the small energetic changes which occur during the rotation of local moments, an 
aspect that makes a proper convergence of the Kohn-Sham iteration very time consuming. Furthermore, our findings 
show that tilted spins due to frustration can even be favored when the total magnetization vanishes (Crg) or the 
corresponding ground state is collinear like in Cr,j. The free variation of the spin quantization axis finally leads to a 
better agreement with the experiment concerning the total magnetization of Cri2 and &13. It can be concluded that 
noncollinear effects appear to be an important ingredient for a deeper understanding of the subtle magnetic properties 
in transition metal clusters. 
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FIG. 1. Geometric and magnetic structures for the energetically lowest noncollinear configurations 
of CrN, 3 < N < 6. The local magnetic moments including their angles with respect to the x- and 
z-axis are indicated by arrows. The interatomic distances are shown in atomic units. 
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M a t 
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M at (refj 


23 \ 
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collinear 


0.0 (-) 


0.0 (-) 
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noncollinear 


0.69 


0.083 


2.0 
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collinear 


0.0 (1.33) 


0.0 (0.12) 


0.0 
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noncollinear 


0.53 


0.054 


0.93 
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noncollinear 


0.0 


0.022 


0.33 
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noncollinear 


0.13 


0.019 


0.29 
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noncollinear 


0.09 


0.015 


0.22 


12 


noncollinear 


0.81 


0.011 


1.67 


13* 


noncollinear 


0.59 


0.008 


1.06 



TABLE I. Total magnetization per atom M a t (in units /ib) for the ground states of 
CrN and gain in binding energy AE nc (in units eV/atom) with respect to their collinear 
counterparts. In case of a collinear ground state the result for the energetically lowest 
noncollinear isomer is given in parentheses. The last column shows the corresponding 
magnetization from the collinear calculation of Cheng and WangH3. The star at N at = 13 
indicates that the structure was not completely relaxed. 
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FIG. 2. 3D-plot of the polar rotation angle 9 (in degrees) from eqn.M) within the plane defined by the 
ionic coordinates of Cr3 (see FigM. We also show the corresponding contour lines with a step size of 10°. 
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FIG. 3. Geometric and magnetic structure of the noncollinear ground state of Cri2. The notation 
follows the one from Fig.pJ. The magnetic moments of the lighter-colored atoms show a dispersion in 
the positive z-direction whereas the darker atoms almost exactly point in the negative z-direction. We 
also display the absolute values of the azimuthal and polar angles of jtax for all 4 groups of atoms. 



